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Abstract 

The reactions 77 — > 7r + 7r~7r + 7r~ and 77 — > 7r + 7r°7r~7r° are studied with the 
L3 detector at LEP in a data sample collected at centre-of-mass energies from 
161 GeV to 209 GeV with a total integrated luminosity of 698 pb _1 . A spin-parity- 
helicity analysis of the p°p° and p + p~ systems for two-photon centre-of-mass energies 
between 1 GeV and 3 GeV shows the dominance of the spin-parity state 2 + with 
helicity 2. The contribution of + and 0~ spin-parity states is also observed, whereas 
contributions of 2~ states and of a state with spin-parity 2 + and zero helicity are 
found to be negligible. 
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1 Introduction 



Several experiments have observed a large cross section near threshold for the reaction 77 — > 
p°p° [1-3]. In contrast, the corresponding cross section for the isospin-related reaction 77 — > 
p + p~ was shown to be small [4,5]. The first spin-parity-helicity analysis of the reaction 77 — > 
7r + 7r~7r + 7r~ was carried out by the TASSO Collaboration [2] by studying angular correlations. 
The data sample consisted of 1722 events for two-photon centre-of-mass energies 1.2 GeV < 
IU 77 < 2.0 GeV. A spin-parity-helicity analysis with higher statistics was performed by the 
ARGUS Collaboration [3] with 5181 events in the region 1.1 GeV < W~ n < 2.3 GeV. Both 
collaborations used similar models and observed the dominance of p°p° states with spin-parity 
J p = 2 + and + . The contribution of negative-parity states was found to be negligible. 

A number of theoretical models [6] were proposed to interpret these experimental results. 
In a t-channel factorization approach [7], the 77 — ■> p°p° cross section is related to photo- 
production and hadronic cross sections at low energies. This model leads to the interpretation 
of the broad enhancement in the 77 — > p°p° cross section around 1.6 GeV as a threshold 
behaviour due to Regge exchange. Other models suggest an s-channel p°p° resonance [8,9], 
either a normal qq state or a four-quark qqqq bound state. In four-quark models, isoscalar and 
isotensor resonances interfere destructively to suppress the 77 — > p + p~ signal and constructively 
to describe the 77 — > p°p° cross section. The proposed models differ substantially in the 
predicted cross section for the production of other vector mesons such as 77 — > p°u and 

77 -> 00. 

This Letter presents the results of a spin-parity-helicity analysis of the reactions 77 — ■> 
7r + 7r~7r + 7r~ and 77 — > 7r + 7r°7r~7r° in data collected by the L3 detector [10] at LEP, using the 
same technique as TASSO and ARGUS. The data samples consist of 7.5 x 10 4 events for the 
e + e~ — > e + e~/T + 7r~7r + 7r~ channel and 7.5 x 10 3 events for the e + e~ — > e + e~7r + 7r°7r~7r° channel. 
These data are selected in the region of quasi-real photons with a maximum virtuality of 
Q 2 ~ 0.02 GeV 2 . The 77 — > p°p° and 77 — > p + p~ cross sections obtained in this analysis are 
compared to the high- virtuality [11, 12] and mid- virtuality [13, 14] data obtained with the same 
detector. 



2 Data and Monte Carlo Samples 

The two-photon production of a p-pair, 77 — > p°p° or 77 — > p + p~, is observed via the reactions 
e + e~ — ► e + e~7r + 7r~7r + 7r~ or e + e~ — ► e + e~7r + 7r°7r~7r°, respectively. Detection of the scattered 
leptons is not required. The data were collected with the L3 detector at e + e~ centre-of-mass 
energies y^s = 161 — 209 GeV, with a total integrated luminosity C e + e - = 697.7 pb _1 and an 
average centre-of-mass energy of 196 GeV. The analysis described in this paper is mainly based 
on the central tracking system and the electromagnetic calorimeter. 

Four-pion Monte Carlo events are generated with the EGPC [15] program. The four- 
momentum of the two-photon system is distributed according to the transverse two-photon 
luminosity function [16]. The pion four-momentum vectors are generated using four-particle 
phase space. The events are then passed through the L3 detector simulation, which uses the 
GEANT [17] and GEISHA [18] programs, and are reconstructed following the same procedure 
as used for the data. 
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3 Event Selection 



The events are collected by two charged-track triggers. The first trigger [19] requires at least 
two wide-angle tracks, back-to-back within ±41° in the plane transverse to the beam. The 
second trigger [20] is based on a neural network which was trained to select low-multiplicity 
events while rejecting beam-gas and beam-wall background. 
Events are selected by requiring: 

• four charged tracks for the e + e~ — > e + e - 7r + 7r~7r + 7r~ reaction and two charged tracks for 
the e + e~ — > e + e~7r + 7r°7r~7r° reaction, with a net charge of zero in each case. A track 
is required to have: more than 12 hits, with at least 60% of possible hits, a transverse 
momentum, p t , greater than 100 MeV and a distance of closest approach to the interaction 
vertex in the transverse plane less than 2 mm. 

• no photons for the 77 — > 7r + 7r~7r + 7r~ reaction and four isolated clusters in the electro- 
magnetic calorimeter for the 77 — ► 7r + 7r°7r~7r° reaction. A photon is defined as an isolated 
shower in the electromagnetic calorimeter consisting of at least two adjacent crystals with 
an energy greater than 100 MeV and with no charged track within 200 mrad. 

• an energy loss dE/dx in the tracking chamber corresponding to the hypothesis that all 
the charged particles are pions, with a confidence level greater than 6%. 

• two pairs of photons each with a good fit to the 7r° decay hypothesis for the 7r + 7r°7r~7r° 
final state. 

To suppress the background from non-exclusive events, the overall transverse momentum 
of the event, |Sp t | 2 , must be less than 0.02 GeV 2 , as shown in Figures ^ and \T]p. The 
resulting samples consist of 74859 and 7535 events for the e + e~ — > e + e~7r + 7r~7r + 7r~ and 
e + e~ — > e + e~7r + 7r 7r _ 7r° reactions, respectively. 

The distributions of the four-pion mass, equal to W 77 for exclusive events, are shown in 
Figures ^ and ^i. The mass resolution is estimated to be 48 MeV for the 7t + 7t~7t + tt~ and 
63 MeV for the 7r + 7r°7r~7r° final states. More than 90% of the events lie in the region 1.0 GeV < 
W 77 < 3.0 GeV, where the spin-parity-helicity analysis is performed. 

The background is dominated by higher-multiplicity final states produced in two-photon 
interactions which are only partially reconstructed. The expected contribution from annihi- 
lation events is negligible. As presented in Figures^ and ^p, the distribution of |Sp t | 2 for 
non-exclusive final states has an exponential form, which is estimated from the data in the 
high |Sp t | 2 region, 0.2 GeV 2 < |Sp t | 2 < 0.8 GeV 2 . Extrapolating this exponential to the signal 
region, |Sp t | 2 < 0.02 GeV 2 , the backgrounds for the and 7r + 7r°7r~7r° final states are 

estimated to be 2.5% and 4%, respectively. 

Figures EH, Efc and EH show the two-dimensional distributions of the masses of 7r + 7r~ combi- 
nations for the selected 7r + 7r~7r + 7T~ events in different W J7 regions. There are two entries per 
event, displayed by ordering the two masses of each entry. Figures Et>,Eli andEF show the tt + ty 
and 7r~7r° mass combinations for the 7r + 7r°7r~7r° channel with two entries per event. 

The two-pion mass resolution is estimated from Monte Carlo simulation to be 25 MeV for 
both the tt + ii~ and ^tt cases. The tt + tt~ and 7r ± 7r° combinations shown in Figure El present 
clear evidence of pp production. For W 71 < 1.6 GeV, the p signal is distorted by threshold 
effects. As increases, the p signal shifts to its nominal mass value, shown by the dotted 
lines in the figure. 
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4 Spin-Parity-Helicity Analysis 

Following the model proposed by the TASSO Collaboration [2], we consider pp production in 
different spin-parity and helicity states (J p , J z ), together with an isotropic production of four 
pions, denoted as "47r" . All states are assumed to be produced incoherently, and therefore no 
interference effects between the final states are taken into account. However, since states of 
different spin-parity and helicity are orthogonal, all interference terms vanish on integrating 
over the angular phase space. Isotropic pirn production, included in previous analyses [3,5], 
corresponds to an unphysical state since C-parity requires the angular momentum between the 
two pions to be odd. We have verified that this state is not essential to reproduce the data. 
The pirn events, if neglected, are absorbed by the An background. 

The analysis is performed in W iy intervals of 100 MeV for 77 — > it + it~it + it~ and 200 MeV 
for 77 — > 7T + 7r 7r~7r°. As pions are bosons, the amplitudes which describe the process must be 
symmetric under interchange of two pions with the same charge and are: 



where m p indicates the mass of the two-pion system and B p (m p ) is the relativistic Breit-Wigner 
amplitude for the p meson [21]. The angular term \1/ jpj z ls{Pii P2) describes the rotational 
properties of the pp state with spin-parity J p and helicity J z . It is constructed by combining 
the spins of the two p mesons, S = si + s 2 , with z projection M s = m\ + m 2 and then adding 
this to the pp orbital angular momentum, L, with z projection M, to obtain the state with 

— * 

total angular momentum J and z projection J z = M s + M: 



where Cf^ t m are the Clebsch-Gordan coefficients, Yz m (£i) are the spherical harmonics and 
£1 = (i?p,Vp), £2 = (1?^+,^+) an d £3 = OV 1 -, ( P 7T +), with -dp and y? p being the polar and 
azimuthal angles of a p meson in the two-photon helicity system. The z axis is chosen parallel 
to the beam direction, which to a good approximation is parallel to the 77 helicity axis. The 
angles ~&^+ and ip n + are the polar and azimuthal angles of the positive pions in the centre-of- 
mass of the first p° meson, with the z axis parallel to the beam axis, the angles 1?^+ and ip n + 
correspond to the second p° meson; for a p~ meson, £ 3 = (•& ir -,(p n -). The indices from 1 to 

4 refer to the four pions using the convention: nfn^n^^ or ti^tt^it^tt^. Since the analysis 
is performed close to threshold, the orbital angular momenta are restricted to L — 0, 1. The 
allowed spin-parity-helicity final states of the pp system in quasi-real two-photon reactions are 
then: ( J p , J z ) = + , 0~, (2+,0), (2 + ,±2) and (2 _ ,0), with the total spin of the pp meson 
system S — 1 or S — 2. States with helicity one are forbidden by helicity conservation and 
spin-one states by the Landau- Yang theorem [22]. 

A maximum-likelihood fit to the data is used in each bin to determine the contributions 
of the four amplitudes: 4tt, + , 0~ and (2 + ,2). The remaining spin-parity states are not 
considered as they have a negligible contribution if included in the fit. 

5 Cross Section 

The cross section for the process k, with fraction determined from the fit, averaged over the 
bin with N events, is 



g.jpj z = B p {m pi )B p (m p2 )^ jp j^sipi, p 2 ) + permutations, 

and 04tt = 1, 




M,mi 
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= 

77 ^ £ e +e- ^fc(W / 77 ) £trg(W" 77 ) / d£ 77 

where J d£ 77 is the two-photon luminosity function integrated over the W 77 bin, e k is the 
selection efficiency and e tTg is the trigger efficiency. The selection efficiencies depend on W 77 as 
well as on the particular wave. They are computed by Monte Carlo simulation, re-weighting 
the events with the amplitudes \gk\ 2 - The efficiencies for the 4-7T process are listed in Tables Q 
and El Similar efficiencies are found for the other processes. The trigger efficiency is studied by 
comparing the response of the two charged-track triggers. The higher-level trigger efficiencies 
are determined using prescaled events. The total trigger efficiency is given in Tables ^ and [21 
The overall efficiencies for the 7r + 7r~7r + 7r~ and 7r + 7r°7r~7r° final states are shown in Figures 
andlHt"- 

The cross sections derived from the fit are presented in Tables^andEl FiguresEfc— f compare 
the total cross sections and the contributions of the individual waves to the 77 — > p°p° and 
77 — > p + p~ processes as a function of the four-pion mass. The 4ir background, which in this 
analysis represents all states which do not correspond to the pp production hypothesis, is similar 
in both channels. It grows from threshold to a value of 20 — 30 nb around 2 GeV and decreases 
toward 3 GeV. In the 7r + 7r~7r + 7r~ channel, the p°p° production has a high cross section, with 
a maximum of about 50 nb at f.6 GeV. It is dominated by the (2 + ,2) state, which has a 
cross section peak-value of about 35 nb. The 7r + 7r 7r~7r° channel, exhibits also a significant 
(2+,2) contribution, but only for 1.6 < IU 77 < 2.2 GeV. Above 1.9 GeV the 77 -> p°p° and 
77 — > p + p~ cross sections are equal within the experimental uncertainties and fall rapidily with 
increasing W 77 . In the + wave, a clear peak is observed in the p°p° channel at W 77 ~ 1.4 GeV, 
perhaps indicative of an s-channel resonance effect, it is absent in the p + p~ channel. The 
77 — > p + p~ cross section peaks near 2 GeV in both + and (2 + , 2) waves, while a similar peak 
is seen for the p°p° in the + wave only. The same qualitative features were observed by the 
ARGUS Collaboration [3,4], which however found a higher peak cross section of ~ 50 nb for the 
(2 + , 2) wave in p°p°. Taking into account the larger experimental uncertainties on the ARGUS 
data, a peak is also seen in the + wave at IU 77 ~ 1.4 GeV. However, at higher mass values, 
IU 77 > 2 GeV, only the much higher statistics of the present experiment are able to provide 
cross section measurements, so no comparison is possible. 

To evaluate the quality of the fit and of the detector modeling we compare several distribu- 
tions of the data with a Monte Carlo simulation normalized to the fit results. Figure 0] shows 
the distributions of p t and of the cosine of the polar angle of the charged or neutral pion closest 
to the beam line. The two-pion mass combinations, tt + tt~ and ^ir , and the production angles 
of the pions in the two-pion centre-of-mass system with respect to the beam direction (Adair 
angle) are plotted in Figure El Four entries per event are considered and the data are plotted 
before acceptance corrections. The agreement with the Monte Carlo simulation is adequate, 
considering the simplicity of the model and the high statistics of the data sample. The angu- 
lar distributions in Figures and 01 are of the general form sin 2 9, indicating a dominantly 
transverse polarisation for the produced p. 

6 Background Estimation and Systematics 

The fraction of non-exclusive background in the pp sample is derived by performing a spin- 
parity-helicity analysis of the background data sample, defined as the region 0.2 GeV 2 < 
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|Xpt| 2 < 0.8 GeV 2 . We find that less than 30% of these events are classified as pp. The 
background contribution in the pp sample is then of the order of 1%. 

Systematic uncertainties on the pp cross sections are due to selection criteria, fitting pro- 
cedures and trigger efficiencies. Uncertainties from the selection procedure are estimated by 
varying the cuts on the quality of the event and on the track definition. They vary between 
3% and 10% for the p°p° channel and between 10% and 24% for the p + p~ channel, depending 
on W 77 . Uncertainties on the model and the fitting procedure are estimated by neglecting in 
turn the + and 0~ waves and including the (2 + ,0) and (2~,0) waves in the fit. Small effects 
from the inclusion of additional spin-parity states are also considered. In the high mass region, 
Wyy > 2 GeV, the contribution of other channels and higher-angular momentum states may 
become important. It was found that such effects could be modelled by including a contribution 
from the isotropic pnir production. In total, these uncertainties for the p°p° and p + p~ channels 
amount to a maximum of 10% and 19%, respectively, for W 77 < 2 GeV and to a maximum of 
60% in the region 2 GeV < W 77 < 3 GeV. Uncertainties on the determination of the trigger 
efficiencies are of a statistical nature and affect mainly the 7r + 7r°7r~7r° channel, where they vary 
between 2% and 6%. They are below 1% for the 7r + 7r~7r + 7r~ channel. Uncertainties on the 
background level are below 1% for both channels. 

7 Discussion 

A spin-parity-helicity analysis of four-pion final states produced in quasi-real two-photon colli- 
sions at LEP benefits from data statistics an order of magnitude higher than previous analyses. 
Several characteristics of the 77 — > p°p° and 77 — > p + p~ processes, which were previously 
observed [3,5], are confirmed: 

• In both channels, the (2 + ,2) wave is dominant. Small but significant + and 0~ waves 
are also observed. 

• The 77 — > p°p° process has a high cross section extending from threshold to about 2 GeV, 
while the cross section of the 77 — > p + p~ process is low in this range. In FigurelBlthe mass 
spectra of the present results are compared to those we obtained at higher Q 2 [11-14]. 
The ratio 

R = Aff ec(p V)/ E A ^c(pV), 

where Acr ee = A£ 77 cr tot (77 — > pp) and the sum is for the region 1.1 GeV < W 77 < 
2.1 GeV, is found to be 

R = 0.42 ± 0.05 ± 0.09 for Q 2 < 0.02 GeV 2 . 

The first uncertainty is statistical and the second systematic, calculated assuming the 
systematic uncertainties for the two processes to be fully uncorrelated. This ratio increases 
with the photon virtuality. At higher Q 2 we previously obtained: 

R = 0.62 ± 0.10 ± 0.09 for 0.2 GeV 2 < Q 2 < 0.85 GeV 2 [14], 
R = 1.81 ±0.47 ±0.22 for 1.2 GeV 2 < Q 2 < 8.5 GeV 2 [14], 
R = 2.2 ±1.1 ±0.6 for 8.8 GeV 2 < Q 2 < 30 GeV 2 [12]. 

These measurements are consistent with the presence of an s-channel enhancement at low 
p°p° mass values which decreases rapidly with Q 2 . If interpreted as an effect of s-channel 
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resonances, the observed ratio between 77 — > p + p~ and 77 — > p p production implies 
the possible existence of an isospin-2 state [8,9]. Such an interpretation of our data was 
recently presented in Reference 23. 

• At higher masses, > 2 GeV, the 77 — > p + p~ and 77 — > p°p° cross sections are equal, 
within the experimental uncertainties. In both cases, the cross section decreases rapidly 
for IU 77 ~ 3 GeV. 

The Q 2 -dependence of the two-photon cross section is presented in Figure [H for the full mass- 
region 1.1 GeV < W yy < 3. GeV. The p°p° cross section exceeds the p + p~ one at low Q 2 while 
a cross-over is observed in the vicinity of Q 2 ~ 1 GeV 2 . A Generalised Vector Dominance fit, 
GVDM [24], which reproduces well all the mid-virtuality and high-virtuality data [14] for the 
77 — » p°p° cross section, lies below the cross section value obtained at (Q 2 ) = 0.001 GeV 2 . A 
p-pole fit, also presented in the Figure, better describes the low-Q 2 region. 
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Table 1: Cross section measurements and fit results for 77 — > 7r + 7r~7r + 7r~ for different W~ n intervals. N is the number of events in 
a bin, f d£ 77 the two-photon luminosity function, e tT g the trigger efficiency and e the selection efficiency. The cross sections for the 
background, An, and for the different spin-helicity waves are given, along with the total 77 — > p°p° cross section. A double dash 
indicates that no significant contribution to the fit is observed. The first uncertainties are statistical, the second systematic. 
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Table 2: Cross section measurement and fit results for 77 — > 7r + 7r°7r~7r° for different W 77 intervals. iV is the number of events in a 
bin, J d£ 77 the two-photon luminosity function, e t r g the trigger efficiency and e the selection efficiency. The cross sections for the 
background, 4n, and for the different spin-helicity waves are given together with the total 77 — > p + p~ cross section. A double dash 
indicates that no significant contribution to the fit is observed. The first uncertainties are statistical, the second systematic. 




Figure 1: Distributions of |Sp t | 2 for a) e + e 
events. The hatched areas represent the estimated non-exclusive backgrounds. The cut values 
are shown by the arrows. Distributions of the four-pion mass for c) e + e~ — > e + e~7r + 7r~7r + 7r~ 
and d) e + e~ — > e + e~7r + 7r°7r~7r° events. Only events within the region indicated by the arrows 
are further analysed. 
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Figure 2: Two-dimensional distributions of two-pion masses in three different W 77 regions. For 
a), c) and e) the 7r + 7r~ combinations from the 7t + 7t~7t + tt~ final-state are shown as low-mass vs. 
high-mass, with two entries per event. In b), d) and f) the 7r + 7r° vs. tt~7t° combinations from 
the 7r + 7r°7r~7T final-state are shown, with two entries per event. The dotted lines indicate the 
nominal mass value of the p meson. 
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Figure 3: 

e + e~7r + 7r 7r _ 7r° processes: c) the total 77 
+ , f) 0~ contributions. The error bars show the statistical uncertainties. 
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Figure 4: Comparison of the Monte Carlo simulation normalised to the fit results to the data: 
a) and b) transverse momentum of the charged or neutral pion closest to the beam line, re- 
specively; c) and d) cosine of the polar angle of the charged or neutral pion closest to the beam 
line, respectively. The statistical uncertainty on the Monte Carlo distributions -not shown- is 
comparable to that of the data. 
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Figure 5: Comparison of the Monte Carlo simulation normalised to the fit results to the data: 
a) two-pion opposite-sign mass combinations for 77 — > 7r + 7r~7r + 7r~ (four entries per event) b) 
two-pion charged mass combinations for 77 — > n + n n^n (four entries per event) c) cos9 w +, 
where n + is the production angle with respect to the beam axis in the unlike-sign two-pion 
centre-of-mass system for 77 — > 7r + 7r~7r + 7r~ (four entries per event) and d) cos^±, where 9 n ± 
corresponds to the 9 n + angle for the 7r + /T°7r~7r° system (four entries per event). The error bars 
show the statistical uncertainties. 
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Figure 6: The 77 — > p°p° and 77 — > p + p cross sections as a function of the four-pion mass 
a) as obtained in the present analysis at Q 2 < 0.02 GeV 2 compared to previous L3 results 
obtained at b) 0.20 GeV 2 < Q 2 < 0.85 GeV 2 [13,14], c) 1.2 GeV 2 < Q 2 < 8.5 GeV 2 [11,12] and 
d) 8.8 GeV 2 < Q 2 < 30 GeV 2 [11, 12]. The error bars show the statistical uncertainties. 
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Figure 7: a) The 77 — > p°p° and 77 — > p + p~ cross sections as a function of Q 2 . b) Comparison 
of the 77 — * p°p° cross section as a function of Q 2 to a GVDM and a simple p-pole form- 
factor dependence, both fitted to previous L3 data at higher Q 2 [14]. The error bars show the 
statistical uncertainties. 
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